Entanglement scaling in critical two-dimensional fermionic and bosonic systems 
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We relate the reduced density matrices of quadratic bosonic and fermionic models to their Green's 
function matrices in a unified way and calculate the scaling of bipartite entanglement of finite systems 
in an infinite universe exactly. For critical fermionic 2D systems at T = 0, two regimes of scaling are 
identified: generically, we find a logarithmic correction to the area law with a prefactor dependence 
on the chemical potential that confirms earlier predictions based on the Widom conjecture. If, 
however, the Fermi surface of the critical system is zero-dimensional, we find an area law with a 
sublogarithmic correction. For a critical bosonic 2D array of coupled oscillators at T = 0, our results 
show that entanglement follows the area law without corrections. 
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Entanglement is a key feature of the non-classical na- 
ture of quantum mechanics. It is a necessary resource 
for quantum computation and is at the heart of inter- 
esting connections between quantum information theory 
and traditional quantum many-body theory, such as in 
quantum critical phenomena 0,12,1^] or the quantum Hall 
effect HHl. 

One of the most widely used entanglement measures 
is the entropy of bipartite entanglement, which is noth- 
ing but the von Neumann entropy of quantum statistics: 
For a pure state I^'.ab) of a bipartite "universe" AB con- 
sisting of system A and environment B it is given by 
Sa = -Trp^logjp^, where pA = Tr^ is 
the reduced density matrix of system A. 

An important question to ask is how entanglement en- 
tropy scales with the size of the system, assuming the 
universe to be in the thermodynamic limit. This was first 
studied by Beckenstein in the context of black hole en- 
tropy As opposed to thermodynamic entropy, which 
is extensive, entanglement entropy was found to be pro- 
portional to the area of the black hole's event horizon, 
its physical locus being essentially the hypersurface sep- 
arating system and environment. Entanglement entropy 
scaling hence depends decisively on the dimension d of 
the universe. 

This observation has given rise to a long string of stud- 
ies of this so-called area law. In one dimension d = 1, 
scaling is well understood both for fermions [2, 7-11] and 
bosons 0, Foi' one-dimensional spin chains, one 

finds that the entanglement entropy Sa{L) of a system A 
of linear size L saturates away from criticality, but scales 
as log2 L at criticality . In the latter case, conformal 
field theory (CFT) yields [ll[ll Sa{L) = ^ log^ L + k, 
where c and c are the holomorphic and the anti-holomor- 
phic central charges of the field theory. Essentially, there 
is no physical limit to the boundary region between sys- 
tem and environment. 

The situation is far less clear in higher dimensions d > 
1. The area law implies that the entanglement away from 
criticality is essentially proportional to the surface area 



of system A 



Sa^L 



(1) 



as confirmed in analytical calculations for non-critical 
bosonic coupled oscillators [l^ . 

At criticality, the correlation length diverges and one 
may expect corrections to the area law, as for d = 1. For 
critical ground states of fermionic tight-binding Hamil- 
tonians entanglement was indeed found to scale as 



Sa^ L''-^\og^L, 



(2) 



for both lattice models ^3 and continuous fields p^ . 
The prefactor could only be derive d 1 181 assuming (i) the 
validity of the Widom conjecture |19l| and (ii) its appli- 
cability to the functional form of binary entropy. For 
bosons at criticality, numerical evidence for the area law 
was found for a three-dimensional array of coupled 
oscillators 20] . Callan and Wilczek derived the area law 
in approximative field theoretical calculations |2 j ]. 

Beyond the fundamental physical interest, entangle- 
ment scaling sets the scope of entanglement-based numer- 
ical methods such as the density-matrix renormalization- 
group (DMRG) [2^ E^ . as the computation time re- 
quired to simulate a quantum state using these methods 
on classical computers increases exponentially with its 
entanglement entropy. 

In this Letter, we study the bipartite entanglement 
in a unified treatment of a class of exactly solvable two- 
dimensional fermionic and bosonic models at T = 0. To 
this purpose, we relate the reduced density matrix of a 
quadratic model to its Green's function matrices, gen- 
eralizing work by Cheong and Henley |2^ based on a 
coherent-state method developed by Chung and Peschel 
(2^. For the critical fermionic two-dimensional tight- 
binding model we find as expected lO, but our exact 
calculation allows to identify the dependence of the scal- 
ing law prefactor on the chemical potential /i. We exactly 
verify the behavior predicted in [iSj, where the validity 
of the Widom conjecture and its applicability to the bi- 
nary entropy were assumed. Interestingly, we observe a 
sublogarithmic correction to the area law if the gap of the 
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model closes in a zero-dimensional region of momentum 
space (i.e. one or more points). For a critical bosonic 
two-dimensional model of coupled harmonic oscillators 
we find the entanglement to saturate to the area law 
which confirms |2Clli21i |. 

The generic quadratic Hamiltonians studied here are 



Hp.B — 



E 



lialBtja] + h.c.) 



(3) 



where = Ci and Ui = bi are fermionic and bosonic 
operators for Hp and Hb respectively. 

Calculating entanglement from Green's matrices. — We 
consider a bipartite universe AB of N levels (or sites). 
System A consists of n sites; in our calculations we will 
eventually take the thermodynamic limit N oo. The 
relation between the Green's function matrices of system 
A and its reduced density matrix = Trg p can be 
derived by determining the matrix elements of the full 
density matrix p with respect to coherent states and in- 
tegrating out the variables of the environment B. Here, 
we focus on the key steps and results of this method. A 
detailed derivation will be given elsewhere. 

(a) Fermionic systems. The block Green's function 



matrix with respect to the operators Ai 



B, 



c] — Ci, as defined by 



[GsAlij = TrpBiAj with i,j e A, 



and 



(4) 



can be obtained exactly for the solvable Hamiltonian Hp, 
Eq. (O, following |2^. It can then be shown that p^ is 
given by 



{^\pa\^')= detail -Gba) 



X e 2 



5(€*-€')^-(GBA + l)(C?i5A-l)-'-(€*+C') 



(5) 



where ^ — {^i • • • f„} are the Grassmann variables as- 
sociated with system A, and |^) are the corresponding 
coherent states with Ci\^) = 

To calculate the (entanglement) entropy of system A, 
we diagonalize p_4 by the Bogoliubov transformation 



fq ^ Si 



(6) 



where PP^ — QQ^ = 1 (due to the anti-commutation 
rules), Pq G%j^ = VqQq and Qq Gba = VqPq ■ The diag- 
onalized reduced density matrix reads 



PA 



(7) 



with pseudo-energies Sq = In jq:^ , such that the entan- 
glement entropy reads 



h{x) 



SA^TUhC-^l 



-xXog^x- (1 -a;)log2(l 



(8) 



(9) 



being the binary entropy. 

(b) Bosonic systems. For the quadratic Hamiltonian 
Hb the block Green's function matrices Gaa and Gbb 
with respect to the operators Ai = b\ + bi and Bi — b\ — bi 
can be obtained as in With respect to the bosonic 
coherent states 6^10) = the reduced density matrix 

then reads 



(10) 



where K' = ^ dci (1 + Gaa){^ — Gbb) is determined by 
the normalization of p^. 

The Bogoliubov transformation 



9q=Y.i 



b] 



(11) 



with P^Q = Q^P = 1, PqGAA = PqQq aud Q qG B B = 

—pqPq diagonalizes p^, giving 



PA 



(12) 



where = In (yjf-^r^ are pseudo-energies. The entan- 
glement entropy Sa is the sum of the quasi-particle en- 
tropies. 



Sa = l°g2 ^ - ^ l0g2 ^ 



(13) 



Critical fermionic entanglement and the Widom con- 
jecture. — The form of the logarithmic correction to the 
entanglement in d > 1 dimensional critical fermion mod- 
els and bounds on it have been derived by Wolf [l7j . 
Gioev and Klich Assuming that the Widom conjec- 
ture ,19] holds also for d> 1 and that the non-analyticity 
of the binary entropy h can be i gno red at one point in 
the calculation, Gioev and Klich 'l8i| arrive at 

Sa = Sn{L) = c{p)L\og^L + o{L\og^L), (14) 



c{p)^^^l dS, [ d5fc|n,-nk|, (15) 
1^ Jan Jdr{p.) 

where ft is the real-space region of A, rescaled by L such 
that Vol(ri) = 1. Vectors nx and nk denote the nor- 
mal vectors on the surface dft and the Fermi surface 
9r(^).With the method introduced above, one can calcu- 
late the entanglement for finite L exactly and thus check 
()14|l . also shedding some light on the validity of the as- 
sumptions leading to H15(l . 

The dispersion relation of the two-dimensional tight- 
binding model with periodic boundary conditions 

H ^ - Y.x,y {<i,yCx+l,v + cl^yC^^y+l + h.c) (16) 
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FIG. 1: The prefactor c{fi) in the entanglement scahng law as 
a function of the chemical potential fj, for the ground state of 
the two-dimensional fermionic tight-binding model in compar- 
ison to the result of Gioev and Klich 1^ . Insets show the hop- 
ping parameters and the Fermi surfaces for /i = — 3,— 2,— 1,0. 



FIG. 2: The scaling prefactor c(/i) for the ground state of 
a two-dimensional fermionic tight-binding model with next- 
nearest neighbor hoppings in comparison to the result of j3- 
Insets show the hopping parameters and the Fermi surfaces 
for fi £ [—0.25, —1.75] in the quartered Brillouin zone. 



is £'(k) = — 2-(cos kx + cos ky). The ground state Green's 
function matrix, from which we calculate the entangle- 
ment, reads in the thermodynamic limit 



Gr 



/l 



r(M) (a^r 



d^k „ik-(r- 



(17) 



with r — {x,y). Fig. ^ shows the scaling prefactor c(/i) 
as fitted from the exact entanglement of an L x L square 
with the rest of the universe, which was obtained from 
(ISJ. It is in excellent agreement with (|15|l and supports 
thus the Widom conjecture for d = 2. The same agree- 
ment was found in the model 

^ = -E.,,((i + (-if)4,,c..,+i 

t t , \ 

+ 4,yCx+l,y+l + C}^,yCx^l.y + l + h.C.) 

which has a two-banded dispersion relation -E(k) — 
±2-y/l -I- 4cosA:2; cos^ ky -f Acos^ cos^ ky and a discon- 
nected Fermi surface for n e [—2, 2], Fig. [21 

Especially for a comparison with bosonic systems, 
it is interesting to investigate models with a zero- 
dimensional Fermi surface. In particular we choose the 
two-dimensional model 



+ {i + {-iy+y)cl^c.,^y+, + h.c.), 



(19) 



which has for < h < 1 the two-band dispersion re- 
lation Eik) = ±2-\/l + cos^ kx + 2h cos k^ cos ky, i.e. 
a gap of size 4(1 — /i) at k = (tt, 0). Fig. |31 shows for 
fi = and /i — > 1 how the entanglement converges to 
the area law with a sublogarithmic correction, Sq{L) = 
L ■ o(log2L), meaning limL->oo Sn{L)/{L\og2 L) = 0. 
The curves Sfi{L)/L for finite gaps were extrapolated 
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FIG. 3: The upper right inset shows the entanglement entropy 
per surface unit Sn{L)/L (block of n = sites) for a two- 
dimensional fermionic tight-binding model with modulated 
vertical hopping (see the second inset) at T = 0. The energy 
gap 4(1 — h) closes in the point k = (tt, 0). The extrapolation 
limL_>oo Sn{L)/L suggests a divergence for /i ^ 1. 



to obtain Yitcil^oo Sn{L) / L. Those values indicate in- 
deed a divergence for ft, ^ 1. This result is consistent 
with Eq. (|14|l . as the scaling coefficient c(/i), Eq. H15|l . 
vanishes for systems in d > 1 dimensions with a zero- 
dimensional Fermi surface. Further investigations have 
to determine the analytical form of the sublogarithmic 
correction and its universality. 

Critical bosonic entanglement. — An important ques- 
tion is whether the logarithmic correction observed in 
the entanglement scaling law for critical one-dimensional 
bosonic systems is also present in higher dimensional sys- 
tems. To investigate this, we examine a two-dimensional 
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system of coupled oscillators 

H = hJ:.^y iKy+ -^o'^lv (20^ 

where ^x,y, ^x.y and loq are coordinate, momen- 
tum and self-frequency of the oscillator at site r = 
{x,y). The masses and coupling strengths are set to 
unity. The system has the dispersion relation E(k) — 

lOq -f 4sin^ fca;/2 -I- 4sin^ ky/2, i.e. a gap at k = 0. 

In the low-energy limit, the harmonic oscillators can be 
reduced to a field theory only containing (Vi/))^, which 
describes a massless free bosonic model. The scaling of 
entanglement in this model has been studied by Srednicki 
|20| numerically in ci = 3 dimensions and by Callan and 
Wilczek |2l| with approximate field theoretical methods 
for all d > 1. Both provide evidence for the area law 

Applying the transformation fej = y^-($i + ^ni) with 
Lo ~ y/ujQ + 4, the Hamiltonian (|2()(l is mapped to the 
canonical form and is thus amenable to the method 
introduced above. The translationally invariant block 
Green's function matrices Gaa and Gbb are for T ~ 

GAA{r, 0) = ^ (Pk -g^ cos kxX cos kyy 

GBB{r,0) = -72 d"^ k^^ cos k^x cos kyy 

and the entanglement is obtained from Eq. (|13|l . Spe- 
cial care has to be taken for the limit loq —> 0, as 
this results in a singularity of the integrand for Gaa- 
Fig. 21 displays the entanglement entropy as a function of 
the linear block size L for several luq. The curves con- 
verge for Wo — > and a finite-size scaling analysis yields 
limijQ^o liniL^oo 'S'q(L)/L ~ 0.45, i.e. the critical model 
obeys for d = 2 the area law Sn{L) « 0.45 • L. 

Conclusions. — A relation between Green's function 
matrices of quadratic fermionic and bosonic Hamiltoni- 
ans to reduced density matrices was used to study bi- 
partite entanglement in critical 2D systems. We iden- 
tified and presented exact quantitative results for three 
different regimes of entanglement scaling. Those find- 
ings demonstrate the subtle nature of entanglement at 
criticality, the physical explanation of which remains a 
challenging topic for future research. 
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